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Abstract 

Background: The left and right amygdalae are key regions distinctly involved in emotion-regulation processes. Individual 
differences, such as personality features, may affect the implicated neurocircuits. The lateralized amygdala affective 
processing linked with the temperament dimension Harm Avoidance (HA) remains poorly understood. Resting state 
functional connectivity imaging (rsFC) may provide more insight into these neuronal processes. 

Methods:\n 56 drug-naive healthy female subjects, we have examined the relationship between the personality dimension 
HA on lateralized amygdala rsFC. 

Results: Across all subjects, left and right amygdalae were connected with distinct regions mainly within the ipsilateral 
hemisphere. Females scoring higher on HA displayed stronger left amygdala rsFC with ventromedial prefrontal cortical 
(vmPFC) regions involved in affective disturbances. In high HA scorers, we also observed stronger right amygdala rsFC with 
the dorsomedial prefrontal cortex (dmPFC), which is implicated in negative affect regulation. 

Conclusions: In healthy females, left and right amygdalae seem implicated in distinct mPFC brain networks related to HA 
and may represent a vulnerability marker for sensitivity to stress and anxiety (disorders). 
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Introduction 

Emotions involve brain networks including (pre)frontal cortical 
and limbic areas [1,2]. Within these emotional networks the 
amygdalae play a crucial role [3,4] . Biologically oriented theories 
suggest specific affective information-processing roles for the left 
and the right amygdala [5,6]. An emotional stimulus automatically 
activates the right amygdala, which is thought to play a role in 
dynamic emotional stimulus detection, while the left amygdala 
seems to be more involved in specific, sustained stimulus 
evaluation [7,8]. However, how individual differences can affect 
left and right amygdala related neurocircuits differently remains 
poorly understood [9-12]. Trait and state anxiety has been found 
to modulate amygdala resting-state functional connectivity (rsFC) 
related to ventromedial prefrontal cortical (vmPFC), but not with 
the dorsomedial prefrontal cortical (dmPFC) activity [13]. 



Together with the amygdalae these brain regions are thought to 
be involved in the neuronal circuits of fear behavior, in self- 
referential processing and social interactions [14,15]. Further- 
more, it has been suggested that in order to stop the generation of 
anxious states the strength of amygdala-mPFC functional 
connectivity during rest represents efficient crosstalk between 
these brain regions [16,17]. 

Only recently, researchers became interested in the relationship 
of amygdala rsFC with personality features [18], such as Harm 
Avoidance (HA). Cloningers' psychobiological theory on person- 
ality and genetic inheritance states that scoring high on the 
temperament factor HA is related to increased behavioral 
inhibition and implies a genetically determined bias towards being 
cautious, apprehensive and overly pessimistic [19]. Healthy 
individuals scoring high on HA are more at risk for developing 
mood- and anxiety disorders in the course of their lives [20,21]. 
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Based on anatomical parcellations of the amygdalae, Li and 
colleagues [22], have reported on sex-related amygdala rsFC 
differences in relation to HA. In spite that functional imaging data 
point to lateralization differences in amygdala emotional func- 
tioning in healthy participants, with especially the left amygdala 
implicated in negative affect [6], and a major topic in our line of 
research [12,23-25], to date it remains unclear whether the 
temperament dimension HA may affect left or right amygdala 
rsFC in relation to the mPFC differendy. Brain imaging 
approaches such as resting-state fMRI combined with HA 
measurements may increase our understanding of how behavioral 
more inhibited individuals with the tendency to be more 
pessimistic could be at higher risk to develop affective disorders 
[17,26]. 

Consequendy, the aim of the current study is to test the 
hypothesis that in a homogeneous sample of females - never 
documented to have suffered from neuropsychiatric illnesses - 
individual scores on HA are related to differential left and right 
amygdala - mPFC coupling. Importandy, the selection of the left 
and right amygdala nodes was based on brain anatomical 
coordinates provided by a neuroimaging study of emotion 
processing and emotion regulation in women resilient or 
susceptible to the depressogenic effects of early life stress [27]. 
These nodes fall within the area referred to as the Superficial 
Amygdala (SA), not surprisingly reported to be involved in the 
processing of social information [28,29] and especially relevant to 
Harm Avoidance. 

Across all subjects, we hypothesized the existence of rsFC 
differences for the connections of left and right amygdala with 
distinct regions in the brain. We hypothesized in high HA scoring 
females stronger rsFC correlations between predominantly the left 
amygdala seed and the mPFC. Within the mPFC, we expected in 
particular left amygdala rsFC-HA correlations with the vmPFC. 
Because amygdala lateralization differences are not consistendy 
reported for the dmPFC, we hypothesized no lateralized amygdala 
rsFC-HA correlations with dmPFC areas. 

Methods and Experimental Procedures 

1. Participants 

The study was approved by the ethics committee of our 
University Hospital (UZBrussel) and in accordance with the 
guidelines laid down in the declaration of Helsinki (2004). All 
participants gave written informed consent. This study was part of 
a larger project investigating several neuro-cognitive markers in 
affective disorders. After the structural MRI, all participants went 
through the rs-fMRI. Hereafter other psychological imaging 
paradigms were performed, not related to the current study. 

Sixty right-handed female individuals (mean age = 21.7 y, 
sd = 2.5), all university students, were recruited. Right-handedness 
was assessed with the van Strien questionnaire [30]. Because 
besides gender also age may confound rsFC results, all participants 
were selected within a narrow age range [31]. Participants taking 
medication, other than birth-control pills, were excluded. None of 
the participants reported to have ever used psychotropic 
medications such as antidepressants, mood stabilizers or antipsy- 
chotics, and all were free of illicit drugs. To exclude psychiatric or 
neurological diseases, all volunteers were screened by the first 
author (C.B). Psychiatric disorders were assessed by the Dutch 
version of the Mini-International Neuropsychiatric Interview 
(MINI) [32]. Participants with a psychiatric disorder and/or a 
score higher than eight on the Beck Depression Inventory (BDI-II 
[33]) were excluded. 



2. Temperament and Character Inventory 

The Temperament and Character Inventory (TCI) is a 240- 
item questionnaire developed by Cloninger and colleagues 
[19,34]. The questionnaire is based on a psychobiological model 
that aims to explain individual differences in personality traits [35]. 
The TCI consists of 4 temperament scales (Harm Avoidance (HA), 
Novelty seeking (NS), Reward dependence (RD), Persistence (P)), 
and three character scales (Cooperativeness (CO), Self-directed- 
ness (SD) and Self Transcendence (ST)) [34] . This inventory has 
been used in a variety of studies examining psychobiological 
substrates of personality, including neurobiological, neuroimaging 
and genetic methods [12,36,37]. We extracted only the temper- 
ament dimension HA for our purposes (minimum score = 0, 
maximum score is 36). 

3. Scanning Procedure 

During the resting state measurements, involving exactly five 
minutes of scanning, all participants were asked to stay awake with 
their eyes closed and to think of nothing in particular. To reduce 
sensory confounds as much as possible, the light in the room was 
dimmed during scanning. After the scan, the participants were 
asked to confirm that they had been awake throughout the scan 
and had complied with the instructions. All resting state fMRI 
scans were performed on Monday afternoons, between 3:00 pm 
and 6:00 pm. 

All scans were performed on a 3T Philips Achieva MRI system 
(Philips, Best, The Netherlands) with an eight channel SENSE 
head coil. fMRI measurement was done using a SE-EPI sequence 
(TR/TE = 3000/70 ms; Flip angle = 90°; FOV = 230x230 mm 2 ; 
resolution = 1.80 x 1.80 mm ; Slice thickness/gap = 4.00/ 
1.00 mm; number of slices = 24; number of dynamics = 100; 
dynamic time resolution = 3000 ms). After the fMRI scan a 3D 
anatomical scan using a 3D Tl TFE sequence (TR/TE = 12.00/ 
3.71ms; Flip angle=10°; FOV = 240 x240 x200 mm 3 ; resolu- 
tion^. 00x1. 00x2. 00 mm' 5 ; number of slices=100) was per- 
formed, yielding an anatomical underlay for the fMRI results. 

The fMRI data were analyzed with the SPM8 software 
(Wellcome Department of Cognitive Neurology, London, UK). 
Slice-time correction was performed to correct for small differ- 
ences in the time offset of consecutively measured slices. Hereafter, 
the images were realigned to the first volume of the time series in 
order to correct for head movements. Subsequently, all fMRI 
brain volumes were normalized to the EPI MNI template; 
resampled to 3-mm isotropic voxels and spatially smoothed using 
an 8-mm full-width half-maximum Gaussian kernel. The anatom- 
ical scans were normalized to the Tl MNI template. 

Several further processing steps preceded the voxel-based 
correlation analysis. Data were linearly detrended and band-pass 
filtered (0.01-0.08 Hz). Spurious or nonspecific sources of 
variance were removed from the data through linear regression 
of: 1) the six head-motion parameters obtained in the realigning 
step, 2) the signal from a region in the cerebrospinal fluid, 3) the 
signal from a region centered in the white matter. As proposed by 
Murphy et al. [38] and Weissenbacher et al. [39] resting state data 
were processed without global signal regression. Correlation maps 
were obtained by extracting the BOLD time course from a seed 
region, then computing the correlation coefficients characterizing 
the correlations between that time course and the time courses 
from all other brain voxels. The seed regions were 6-mm-diameter 
spheres designed to encompass the left (MNI coordinates x= —20, 
y=— 4, z=— 15) or right amygdala (x = 22, y= — 2, z=— 15). 
These MNI coordinates were selected following the recent paper 
of Cisler et al [27]. To combine results across subjects and 
compute statistical significance, Fisher's r-to-Z transformation was 
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Figure 1. Left and right amygdala rsFCs. Transversal slides displaying the results of A) the one-sample t-test for the left-amygdala rsFC and B) 
the one-sample f-test for the right-amygdala rsFC. Colors from yellow to red represent significantly stronger FC and colors from green to blue 
represent the opposite. C) the comparison of the rsFCs of left vs. right amygdala the paired f-test between the rsFCs of the left and right amygdala 
seeds. Colors from yellow to red represent significantly stronger FC with the left than with the right amygdala. Colors from green to blue represent 
the opposite: significantly stronger FC with the right amygdala than with the left. The amygdala seeds are displayed by a white circle. For an overview 
of all significant clusters see Table 1, P = posterior, L = left, R = right. 
doi:1 0.1 371 /journal.pone.0095740.g001 




Figure 2. Amygdala rsFC -HA correlation analysis. Transversal slide exhibiting positive (yellow to red) and negative (green to blue) correlation 
clusters for the correlation between the rsFC of the left amygdala seed (Left crosshair on white sphere; MNI coordinates: x= —20, y = -4, z = -15), the 
right amygdala seed (Right crosshair on white sphere; MNI coordinates: x = 22, y = - 2, z = -15), and HA. For an overview of all significant clusters see 
Table 1 and Table 2. P = posterior, L = left, R = right, BA = Brodmann area. 
doi:1 0.1 371 /journal.pone.0095740.g002 
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used to convert these correlation maps into Z maps (maps 
quantifying local 'rsFC strength', or simply 'rsFC'). The Z maps 
were submitted to a random-effects analysis in SPM8. A one- 
sample f-test containing age as covariate was performed for the left 
and right amygdala rsFC separately. To evaluate significant 
differences between left and right amygdala rsFC, a paired i-test 
was performed with age as covariate. All analyses used a cluster 
significance level of p<0.05, corrected for multiple comparisons 
(Family Wise Error (FWE)). We listed all significant clusters with a 
cluster extent threshold (K) of at least 50 voxels. 

Concerning the influence of the temperament dimension HA on 
left and right amygdala rsFC separately, we calculated Pearson's 
correlation coefficients between the Fishery-transformed rsFC 
strength and HA scores for each voxel, producing another set of r- 
maps. To examine our primary research question; the influence of 
HA on lateralized amygdala rsFC, we calculated the Pearson's 
correlation coefficient between the difference of the Fisher-z- 
transformed rsFC strength and HA scores for each voxel left vs. 
m right amygdala rsFC. After Fisher-z transformation on r-maps, we 

S mapped the voxels with/>-values<0.05. The anatomical labels and 

™ Montreal Neurological Institute (MNI) coordinates were obtained 

| by the xjView MATLAB toolbox (http://www.alivelearn.net/ 

"§ xjview). 
m 

% Results 

c 
o 

£ The range in HA scores was between 2 and 30 (mean HA 

0 score = 15.52, sd= 7.04). The Shapiro- Wilk normality test showed 
that HA scores were normally distributed (p = .46). No volunteer 

-g stated to have fallen asleep during scanning. Due to exceeding 

™ 1.5 mm and 1.5 degree in maximum head motion, four female 

volunteers were removed from rs-fMRI analyses, leaving a total of 

56 participants. 

| 1. Amygdala rsFC 

"£ 1.1. Left amygdala rsFC. See Fig. 1 A. The result of the 

| one-sample /-test for the left-amygdala rsFC showed one large 

■| significant cluster in the left parahippocampal gyrus (K = 6285; 

1 MNI coordinates: x=-18, y=-3, z=-15). On the left 
i= hemisphere, this rsFC region included hippocampus, insula and 
o subgenual anterior cingulate cortex (sgACC), putamen and 
.1 claustrum, fusiform gyrus and culmen. This cluster also extended 

0 to the right hippocampus, the left and right gyrus rectus, and the 
<d thalamus. 

Z In addition, the one-sample /-test showed a significant inverse 

£ correlation between the left amygdala and the right middle 

1 occipital gyrus (BA 18: K = 392; MNI coordinates: x=15, y = 

1 -87, z = 6). 

v 1.2. Right amygdala rsFC. The result of the one-sample t- 

2 test for the right-amygdala rsFC showed two significant clusters, 
c One large cluster was situated in the right parahippocampal gyrus 
2 (K = 3376; MNI coordinates: x = 18, y = - 3, z = - 15). A second 
g 5 cluster was located in the right anterior cingulate cortex (BA 24; 
^ | K= 148; x = 3, y = 27, z= 18). See also Fig. 1 B. 

*- ft The one-sample Hest revealed an inverse correlation between 

2 g the right amygdala and the right fusiform gyrus (K = 89; x = 5 1 , 

g- £ y = - 18, z = -30), the left insula (BA 13; K= 62; x = -42, y = -6, 

2 5 z = 0), the left middle frontal gyrus (BA 9; K = 410; x= -33, 

t 1 y = 33, z = 42), and the right (K = 222; x = 48, y = - 66, z = 6) and 

sg ^. left middle temporal gyrus (K = 186; x = —45, y = —69, z = 9). 

^ £ 1.3. Comparison between the left and right amygdala 

ra 5 rsFCs. See also Fig. 1 C. The paired /-test revealed that the 

o is contrast (left amygdala rsFOright amygdala rsFC) yielded a 
significant cluster in the left parahippocampal gyrus, with the 
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Figure 3. Left vs right amygdala rsFC -HA correlation analysis. Left column) Transversal slides displaying the results of the rsFC-HA correlation 
analyses for the left vs. right amygdala rsFC seeds. Colors from yellow to red represent significantly stronger FC with the left compared to the right 
amygdala. Colors from green to blue represent the opposite: significantly stronger FC for the right compared to the left amygdala. For an overview of 
all significant clusters see Table 3. P = posterior, L = left, R = right, BA = Brodmann area. Right column) Scatter plots representing left vs. right amygdala 
rsFC-HA correlations with their respective correlation coefficients. The red circles represent more left vs. right amygdala rsFC; the blue circles 
represent the reverse: more right vs. left amygdala rsFC. 
doi:1 0.1 371 /journal.pone.0095740.g003 



maximum peak in the left amygdala (K = 2889; MNI coordinates: 
x= —21, y= —3, z = —15). Other peak areas were located at the 
left side: in the hippocampus, insula (BA 13), putamen, fusiform 
gyrus, and pons. 

The paired i-test for the contrast (right amygdala rsFOleft 
amygdala rsFC) revealed a significant cluster one in the right 
parahippocampal gyrus with the maximum peak in the right 
amygdala (K = 578; x = 21, y=— 3, z = — 15), hippocampus, 
culmen, and lingual gyrus (BA 17). 

2. Amygdala rsFC-HA correlation analyses 

2.1. Left amygdala. The results of the rsFC-HA correlation 
analysis for the left amygdala seed (MNI coordinates x=— 20, 
y = — 4, z = — 15) yielded a positive association in the left occipital 
gyrus and the sgACC (BA 25), the right amygdala and larger parts 
of the cerebellar regions. Negative correlations were observed in 
the left cerebellum and the right superior frontal gyrus (BA 9). See 
Table 1 and Fig. 2. 

2.2. Right amygdala. The rsFC-HA correlation analysis for 
the right amygdala seed (MNI coordinates: x = 22, y= — 2, z = 



— 1 5) showed a positive association with the left inferior frontal 
gyrus (BA 10) and the post cinglulate gyrus (BA 31). Strong 
correlations were also observed in the right middle frontal gyrus 
(BA 9 as well as BA 11), the right parietal lobe (BA 7), and 
cerebellar regions. Bilateral positive correlations were observed for 
the parietal and occipital cortex and thalamus. Negative correla- 
tions were found in the left ACC (BA 24), the right parahippo- 
campus and cerebellum bilaterally. See Table 2 and Fig. 2. 

2.3. Left vs. right amygdala confined the entire 
sample. rsFC-HA correlation analyses revealed that the high 
scorers on HA displayed stronger left compared to right amygdala 
rsFC within some clusters situated around the left premotor cortex 
(BA 6) and supplementary motor area (SMA). Further, healthy 
females scoring higher on HA displayed stronger left vs. right 
amygdala rsFC in the right inferior frontal gyrus and the right 
sgACC (BA 25). See Table 3 and Fig. 3. 

On the other hand, when comparing right vs. left amygdala 
rsFC, higher scores on HA showed stronger rsFC-HA correlation 
within the right parahippocampal gyrus and prefrontal cortex, 
including the middle frontal gyrus (BA 9), insula, and orbitofrontal 
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Table 3. Results for the correlation between the individual scores on Harm Avoidance and the rsFC of the left vs. right amygdala 
seed. 




Seed Hemisphere 


Cluster size 


Anatomical region 


BA 


Z-value 


Peak coordinates 
(x,y,z) (mm) 


LefORight 

Amygdala 

Correlation 


Left 


81 


Cerebellum posterior 




0.40 


-12 -81 -21 




55 


Precentral gyrus 


6 


0.52 


-18 -12 48 




35 


Precentral gyrus 


6 


0.45 


-12 -21 75 




21 


Supplementary motor area 


6 


0.42 


-15 -3 66 




22 


Middle frontal gyrus 


10 


0.34 


-39 60 -6 


Right 


63 


Inferior frontal gyrus 
(Parahippocampus) 


47 


0.49 


27 9 -24 




21 


Anterior cingulate 


25 


0.37 


3 3-6 




21 


Inferior frontal gyrus 


45 


0.36 


60 24 24 


Right>Left 

Amygdala 

Correlation 


Left 


47 


Posterior cingulate 


30 


0.38 


-3 -54 12 




30 


Medial frontal gyrus 


10 


0.40 


-9 51 15 


Right 


228 


Middle frontal gyrus 


9 


0.52 


39 39 30 




162 


Precentral gyrus 




0.57 


33 -21 42 




49 


Inferior frontal gyrus 


Insula 


0.36 


42 24 -6 




29 


Fronterior superior orbital gyrus 


11 


0.39 


24 42 - 1 5 




25 


Parahippocampal gyrus 




0.35 


18 -12 -24 





For each cluster, we reported the Z-value and MNI coordinates at the position of the maximum, the cluster size (K) and the corresponding Brodmann area (BA). 
doi:1 0.1 371 /joumal.pone.0095740.t003 



cortex (BA 1 1). In this contrast stronger left hemispheric rsFC-HA 
correlations were found in left posterior cingulate (BA 30) and 
medial prefrontal gyri (BA 10). 

2.4. Left vs. right amygdala confined to high HA 
scorers. To evaluate possible involvement of HA in the 
development of affective disorders, we selected those females 
scoring high on HA according Dutch and Flemish normative data 
set (n= 1041) (The Netherlands are a neighboring country closely 
related to Flanders, Belgium). According this data set, which 
provide normative TCI data for males and females separately [40], 
from our 56 female participants 18 scored high or very high on the 
temperament dimension HA. For an overview see Table 4 and 
Fig. 4. 

In short, left vs. right amygdala rsFC showed stronger rsFC-HA 
correlations with the right hemisphere, such as the right insula, but 
importandy also with a larger cluster within the vmPFC, 
comprising both amygdalae and the sgACC, the inferior frontal 
and rectal gyrus, extending to both parahippocampi. This contrast 
also revealed significant rsFC-HA correlation with the dmPFC, 
more in particular the medial prefrontal gyrus (BA 32/24). 

On the other hand, the right vs. left amygdala rsFC showed 
stronger rsFC-HA correlations with the left hemisphere, including 
the basal ganglia, the left superior frontal (BA 9) and bilateral 
medial frontal (BA 10) gyri. 

Discussion 

Although not the main scope of the current research, our overall 
amygdala rsFC observations without the inclusion of the HA 
scores are in line with the findings of Roy and colleagues [41] 



where spontaneous activities in the amygdalae predicted sponta- 
neous activity in similar parahippocampal and prefrontal regions, 
the thalamus, and occipital cortex. Our general rsFC results point 
to distinct functional network connections largely within the same 
hemisphere for left or right amygdala seed separately. Further, the 
amygdala rsFC-HA correlations showed mosdy positive associa- 
tions with temporal, parietal, occipital and cerebellar cortices. 
These areas play critical roles in the perceptual processing of 
socially and emotionally relevant visual information, especially in 
non-clinical samples with higher trait anxiety [42-46]. 

As hypothesized, left vs. right amygdala rsFC-HA correlation 
analyses showed that females scoring higher on HA displayed 
stronger left amygdala FC within mPFC regions. Within the 
vmPFC, more in particular the sgACC, this area is related to 
arousal processes and implicated in a corticolimbic neurocircuit 
associated with 'visceromotor' functions playing an important role 
in modulating affect, such as sadness activation and ruminative 
thought patterns [47] . This functional amygdala - sgACC coupling 
has also been reported in female anxiety patients scoring higher 
than controls on HA [48]. However, our findings seem to be in 
disagreement with the study of Kim and colleagues [13] where 
reverse amygdala -vmPFC FC results were reported in relation to 
higher scores on the State and Trait Anxiety Inventory self-report 
questionnaires (STAI-S, STAI-T [49]). Of note, albeit higher HA 
scorers may display higher anxiety levels, HA and STAI scales do 
not measure the same construct [20] (Cloninger et al., 2006). 
Further, although Kim et al [13] defined the sgACC as part of the 
vmPFC, on a functional level it may be that not the sgACC but the 
more ventral-rostral portions of the ACC and vmPFC are involved 
in regulating strong emotional responses [50-52]. Indeed, in the 
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Figure 4. Left vs right amygdala rsFC -HA correlation analysis confined to high HA scoring females. Transversal slides displaying the 
results of the rsFC-HA correlation analyses for the left vs. right amygdala rsFC seeds. Colors from yellow to red represent significantly stronger FC with 
the left compared to the right amygdala. Colors from green to blue represent the opposite: significantly stronger FC for the right compared to the left 
amygdala. For an overview of all significant clusters see Table 4. 
doi:1 0.1 371 /journal.pone.0095740.g004 



selected group of High HA scorers in particular the left amygdala 
and the right medial frontal gyrus (BA 32/24) were functionally 
connected, the latter pregenual anterior cingulate cortex (pACC) 
documented to control emotional neuronal processes, such as self- 
conscious emotion [53,54]. 

The left vs. right amygdala rsFC-HA correlation analyses 
showed positive right amygdala rsFC-HA with predominantly the 
right hippocampus and insula, the dorsolateral prefrontal cortex 
(DLPFC) and orbitofrontal cortex (OFC). The latter are part of the 
more dorsal parts of the mPFC [17]. See also Table 3. In spite that 
we did not hypothesize lateralized amygdalae rsFC-HA correla- 
tions with dmPFC areas, these findings concur with a right 



prefrontal involvement in the regulation of negative affect. Besides 
that the DLPFC and the amygdala are indirecdy implicated in 
top-down/bottom-up emotion-regulation processes [55], the right 
DLPFC in particular seems to be implied in behavioral inhibition, 
negative affect regulation, increased vigilance and sustained 
attention, uncertainty and ambiguity [56] (Shackman et al., 
2009). Of interest, the right OFC shows co-activations with insular 
parts associated with interoception and gustation [57]. Interocep- 
tive information such as visceral sensations implicated in processes 
of awareness and experiences of aversive responses are thought to 
be channeled into in the right anterior insula [58-60] . Indeed, the 
neurobiological modulation of stress responses has been reported 
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Table 4. Results for the correlation between the selected individual scores high 
left vs. right amygdala seed. 


on Harm Avoidance (n = 


18) and the rsFC of the 




Seed Hemisphere 


Cluster size 


Anatomical region 


BA 


Z-value 


Peak coordinates 
(x,y,z) (mm) 


Left>Right Amygdala 
Correlation 


Left 


286 


Limbic lobe 


- 


0.08 


0 0 -24 




47 


Cerebellum 


- 


2.17 


-51 -63 -24 




30 


Inferior temporal gyrus 




0.98 


-66 -45 -18 


Right 


195 


Medial frontal gyrus 


6 


1.56 


48 0 51 




185 


Insula 


13 


-1.35 


36 -3 9 




181 


Medial frontal gyrus 


32/24 


1.26 


6 15 45 




127 


Insula 


13 


0.77 


27 33 15 




90 


Occipital lobe 


18 


-4.41 


9 -96 12 




79 


Parietal lobe 


40 


-2.83 


63 -48 24 




63 


Occipital lobe 


18 


-3.09 


24 -90 -15 




37 


Middle frontal gyrus 


11 


-1.28 


36 39 -12 




36 


Middle frontal gyrus 


8 


2.00 


36 24 45 




33 


Cerebellum 




1.09 


12 -30 -21 




31 


Inferior frontal gyrus 


45 


-2.94 


60 24 21 




26 


Lentiform nucleus 




-0.60 


18 0 -3 




20 


Cerebellum 




1.09 


30 -66 -42 


Ftight>Left Amygdala 
Correlation 


Left 


252 


Caudate nucleus 




1.32 


-21 -15 24 




97 


Cerebellum 




3.32 


-6 -51 -12 




92 


Superior frontal gyrus 


9 


-0.28 


-18 42 27 




89 


Medial frontal gyrus 


10/11 


-0.003 


-3 54 -9 




76 


Cingulate gyrus 


31 


-0.30 


-12 -27 39 




44 


Precentral gyrus 


6 


1.22 


-51 -9 45 




30 


Precuneus 




3.46 


- 1 2 -42 45 




28 


Thalamus 




-1.51 


-3 -27 -9 




25 


Parietal lobe 


7 


1.05 


-24 -60 48 




21 


Culmen 




3.36 


-33 -33 -30 




20 


Inferior temporal gyrus 




2.31 


-33 3 -36 


Right 


66 


Superior frontal gyrus 


10 


-2.34 


24 51 27 




24 


Supplemental motor area 


6 


0.39 


3 -12 78 





For each cluster, we reported the Z-value and MNI coordinates at the position of the maximum, the cluster size (K) and the corresponding Brodmann area (BA). 
doi:1 0.1 371 /joumal.pone.0095740.t004 



to be lateralized to the right prefrontal cortex [61,62]. Interest- 
ingly, in the selected group of high HA scorers, also the left 
amygdala showed significant FC with the right insular regions. 
Furthermore, the right DLPFC, OFC and insula were found to be 
activated in anticipation to withdrawal-related emotional experi- 
ences [63]. Our current results further provide insight that in 
behavioral more inhibited and cautious individuals, not only the 
right amygdala may play a key role in regulating these processes, 
but in more stress sensitive individuals both amygdalae seem to be 
involved. In addition, rsFC-HA correlations showed that the both 
amygdalae were significantly stronger functionally correlated with 
the medial frontal gyrus (BA 10) also part of the dmPFC area. 
These rostral parts of the dmPFC are associated with emotion 
regulation, sustained attention, memory, and mentalizing process- 
es [64-67]. Being part of the DMN, the BA 10 is implicated when 
individuals make self-relevant affective decisions [68,69]. Impor- 



tantly, right amygdala also correlated with the posterior cingulate 
gyrus, part of the DMN as well, and together with the mPFC and 
hippocampus is thought to be implicated in the processing of 
autobiographical memory, past self-relevant stimuli and future 
prospection [70]. As this DMN is especially implicated when at 
risk for clinical depression, it is tempting to speculate that this right 
amygdala rsFC-HA dmPFC association may represent a 'neuronal 
network vulnerability' for the development of mood disorders in a 
later stage of life. 

Finally, for the higher HA scorers bilateral amygdala rsFC-HA 
correlations extended from the vmPFC to the basal ganglia. This is 
an important observation because the involvement of dopaminer- 
gic nuclei is not surprising. Besides that the amygdala, hippocam- 
pus, and these ventromedial prefrontal cortical areas are key brain 
regions that not only modulate emotions and cognition but also 
the response to stress itself - resulting in hypertrophy of dendritic 
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arborization and increases in spine density [71-73] - the 
mentioned vmPFC areas are consistendy involved in positive 
and negative reward processing (for an overview see Liu et al. 
[74]). This is of particular importance in more behavioral inhibited 
and pessimistic individuals. These dopaminergic neurons coming 
from the ventral tegmental area (VTA) are crucial for the 
recognition of rewards and their consumption [75]. Again, as 
the selected amygdala nodes fall within the area referred to as the 
Superficial Amygdala (SA) which has been shown to be specialized 
in the processing of social information [28-29] , our results add to 
the assumption that individuals scoring high on HA not only 
display more behavioral inhibition and pessimism, but may also be 
more vulnerable to stressful interpersonal experiences. 

Although the selection of psychopathology-free female subjects 
can be considered a major advantage of the study, including only 
healthy women within a certain age range means that we cannot 
generalize our findings to other populations. Because no cardiac 
and respiratory data were collected during rs-fMRI, this should be 
noted as a limitation of our study. As it has been reported that the 
different subnuclei of the amygdalae may have specific functional 
connections with distinct parts of the brain [41,76], by not 
examining dedicated seeds in these subnuclei, important informa- 
tion could have been missed. However, our main research 
objective was to examine rsFC differences in relation to specific 
left and right-sided amygdalar nodes which were documented to 
be involved in emotion regulation brain networks among 
individuals resilient or susceptible to the depressogenic effects of 
early life stress [27]. This makes the choice of these selected nodes 
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particularly relevant in relation to personality features such as 
harm avoidance. And again these nodes fall within the area 
referred to as the Superficial Amygdala involved in social 
information processing [28—29]. Nevertheless, future research 
examining amygdala rsFC in relation to personality features may 
do well to include a larger number of seeds, comprising the 
different amygdalar subnuclei. 

In conclusion, amygdala rsFC analyses in relation to individual 
differences in HA may prove to be a valid method to investigate 
behavioural inhibition and pessimism, possible risk factors for 
mental illness development. Our rsFC-HA results confirm the 
right amygdala's key role in right anterior hemisphere cross-talk in 
females who are likely more stress-sensitive. Furthermore, the 
combination of enhanced left amygdala -vmPFC and right 
amygdala-dmPFC coupling may represent a vulnerability marker 
for females with an elevated risk to develop mood and anxiety 
disorders. Longitudinal follow-up studies in both genders are 
needed to substantiate such hypotheses and to demonstrate 
whether or not such amygdala rsFC-HA patterns within the 
medial prefrontal cortex are of use to predict the development of 
mood and anxiety disorders. 
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